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ABSTRACT 

EX Lup is the prototype of the EXor class of eruptive young stars. These objects 
show optical outbursts which are thought to be related to runaway accretion onto the 
star. In a previous study we observed in-situ crystal formation in the disk of EX Lup 
during its latest outburst in 2008, making the object an ideal laboratory to investi- 
gate circumstellar crystal formation and transport. This outburst was monitored by a 
campaign of ground-based and Spitzer Space Telescope observations. Here we modeled 
the spectral energy distribution of EX Lup in the outburst from optical to millimeter 
wavelengths with a 2D radiative transfer code. Our results showed that the shape of the 
SED at optical wavelengths was more consistent with a single temperature blackbody 
than a temperature distribution. We also found that this single temperature compo- 
nent emitted 80-100 % of the total accretion luminosity. We concluded that a thermal 
instability, the most widely accepted model of EXor outbursts, was likely not the trig- 
gering mechanism of the 2008 outburst of EX Lup. Our mid-infrared Spitzer spectra 
revealed that the strength of all crystalline bands between 8 and 30 fim increased right 
after the end of the outburst. Six months later, however, the crystallinity in the 10 fim 
silicate feature complex decreased. Our modeling of the mid-infrared spectral evolution 
of EX Lup showed that, although vertical mixing should be stronger during the outburst 
than in the quiescent phase, fast radial transport of crystals (e.g., by stellar/disk wind) 
was required to reproduce the observed mid-infrared spectra. 
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1. Introduction 

The EXor class of young eruptive stars, named after tiie prototype EX Lup, was established 
on the basis of the peculiar optical behaviour of these sources (|Herbig] 119771 ) . EXors are known to 
show repetitive optical brightenings separated by faint quiescent periods. Typical time-scales of the 
bright states range from several months to few years. The time between two subsequent outbursts 
is typically several years. During the outbursts the brightness of the source increases by 2-5 mag in 
the V-band. In the quiescent phase the optical spectra of these stars resemble that of typical low- 
mass T Tauri stars with a spectral type of late K to early M. In the bright state, howe ver, the optica l 
spectra are dominated by broad emission lines usually seen in very early- type stars (lHerbigll2008l ). 
The excess emission above the stellar photosphere at infrared wave lengths is also a characteris tics 
of the class which is attributed to a protoplanetary accretion disk ([Hartmann &: KenvonlllQQq ). 



In general, optical o utbursts of EXors are thoug ht to be related to runaway accretion from the 
disk onto the young star (jHartmann &: Kenyonlll996l ) . This scenario is supported by observations of 
accretion tracers (e.g. spectral lines, veiling) which indicate higher m ass accretion rate in the high 
phase compared to the quiescent periods (e.g. iLorenzetti et al.l 120071 ). Little is known, however, on 
the details of the mechanism(s) leading to the increase of the accretion rate. It is generally believed 
that EXor outbursts are caused by the same kind of disk instabilities as those drivin g the eruptions 
in FUOri type objects (EXors are sometimes also called sub-FUors) (jHerbigi 120081 ). 



FUor and EXor outbursts can also play an important role in the thermal processing of proto- 
planetary dust. During outbursts, the enhanced accretion and irradiation luminosity in the central 
regions of the system should lead to a strong increase of the dust temperature in the disk. Wherever 
the temperature exceeds a threshold value, the initially amorp hous silicate grains cr ystallize rapidly. 



Crystalline silicates are not present in the interstellar medium (jKemper et al.ll2005l ) , from where the 



dust content of protoplanetary disks originates. Sharp features of silicate crystals at mid- infrared 



Meeus et al. ( 


2009) 


(e.g., 


Wooden et al. 



and references therein). Although these observations provide evide nce for crystal formati on in pro- 



toplanetary disks, the process itself was not actually observed until [Abraham et al.l ( 20091 ) reported 



the first observation of in-situ crystal formation in its recent outburst of EXLup. 

Crystalline silicates are also ob served in the cold (T<300 K) outer disk, where the temperatur e 
is too low for crystallization (e.g., iBouwman et al.l l2008l : IWatson et al.l l2009l : iMeeus et al.ll2009l ). 
Cometary bodies, which can contain high abundance of crystalline silicates, are also formed in this 
cold part of protoplane tary disks. An im portant question is whether crystals were formed in situ 



in the outer disk (e.g. iDesch et al.l l2005l ) or were formed in th e central r egion of the system at 



high temperature and subsequently transported outwards (e.g. iGaill l2004l ). In the lack of time- 
resolved observations it is hard to distinguish between the two scenarios. In the case of EXLup, 
however, we know that the crystals are formed in the hot (T> 1000 K) inner regions of the system 
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( Abraham et al.ll2009l ). Thus, the recent outburst of EXLup provides the unique opportunity to 
use the sihcate crystals to trace the transport mechanisms in the disk of EX Lup. 



EXLup is a young (3Myr) low-mass (0.6 Mq) MO star (jGras- Velazquez h Ravll2005l ) located 
in the Lupus cloud complex. In the 1950s, EXLup had a large outburst [/SV ~5mag), but only 
smaller amplitude brightenings, separated by few years of quiescent periods, were observed ever 
since. EX Lup has an infrared excess above the stellar photosphere longwards of about 2 /xm, which 
originates in a pro toplanetary disk. So far no evidence for an envelope has been found and indeed 
Sipos et al.l ((20091) successfully modeled the quiescent SED of EX Lup using a protoplanetary disk 
model only. Apart from a large inner hole in the dust disk, no other observational difference has 
been found which would distinguish EX Lup from other T Tauri stars. The 2008 outburst of EX Lup 
was observed with various ground-based instruments as well as with the Spitzer Space Telescope. In 
this paper we analyze these data in order to constrain dust processing and the outburst mechanism. 



2. Observations 

In our observation campaign we used a variety of ground- and space-based telescopes/instruments. 
We attempted to collect as co-temporary measurements as possible to build the SED of EXLup 
in the outburst. In the following we will describe the details of the data reduction for each used 
instrument, while the log of our observations as well as the resulting flux densities are listed in 
Tab.E 



2.1. Optical light curve 

The optical brightness variations during the course of the outburst were monitored by amateur 
and professional astronomers. We collected visual brightness estimations from Albert Joiies, th e 
AAVSOEI database and V-band CCD observations from the ASAS-3 database (jPojmanskil l2002l ). 
For the visual magnitude estimation we used an uncertainty of 0.3 mag, which was the average 
difference between EXLup and the comparison stars. From the ASAS-3 database we used only the 
highest quality measurements (quality flag 'A'), for which the uncertainty of the measured V-band 
magnitude was less than 0.1 mag. The optical light curve is bpresented in Fig.[TJ 



2.2. GROND 

We observed EX Lup with the Gamma- Ray Burst Optical and Near Infrared Detector (GROND), 
which is a 7-channel simultaneous optical-infrared imager at the ESO 2.2 m telescope, on 20 April 
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2008. Optical images, covering a field of view of 5.4' x 5.4', were obtained in the g' ,r',i' ,z' bands. 
The total observing sequenced consisted of 4 dithered exposures with 1 s integration time. At near- 
infrared wavelengths (JHK5) images were obtained in 4 (JH) or 24 (Kg) dither positions, with a 
total of 40 sec exposure time. The field-of view was 10' x 10'. The optical images were reduced 
with standard IRAF tasks, and the final magnitudes were det ermined with apert ure photometry 
and adopting the nominal zeropoint of each optical band from (jGreiner et al.ll2008l ). Near-IR mag- 
nitudes were determined by averaging the photometric zeropoints derived for several stars in the 
field using the 2MASS catalog. Because the near-infrared images were defocused, a large aperture 
was adopted in order to include all flux. The estimated photometric accuracy is 2.5, 2.5 and 5% 
in the J, H, and band, respectively, while in the optical the uncertainty is significantly larger 
(15 %) due to the lack of dedicated calibration measurements. 



2.3. WFI 

EXLup was observed with the Wide Field Imager (WFI) on the 2.2m telescope in La Silla 
on 20 April 2008, with about one hour difference with respect to the GROND observations. The 
field of view of WFI spans 33'x34' with small holes between the 8 CCD detectors. EXLup was 
observed using standard Johnson UBV filters. A total of 3 exposures were taken with each filter, 
with individual exposure times of 4 s for U and 0.5 s for BV. The data was reduced using standard 
IRAF procedures. 

The photometry was done as standard aperture photometry with a 15 pixel aperture. The 
calibration and aperture corrections were obtained from observations of the Landolt fields SA107 
and SAllO (Landolt 1992) observed at 6 different airmasses, resulting in zeropoint errors <0.01 mag. 

Optical photometries corrected for emission line contamination are also listed in Tab.[TJ Cor- 
rection was done by convolving optical spectra (Sicilia-Aguilar et al., in prep) with the WFI and the 
GROND filter curves and calculating the ratio of the total observed emission over the continuum. 
The emission line corrected photometries were used during the SED fitting, since our model did 
not include optical emission lines. 



2.4. NTT SOFI 

EX Lup was observed on 20 April 2008 using SOFI at the ESO 3.5 m New Technology Telescope 
(NTT). We obtained near infrared images using the narrow band filters NB1.215, NB1.71 and 
NB2.195. The observations were performed using the template img-obs-Auto Jitter, taking five 
offsets, and repeating the exposure several times at each jitter position. The high brightness of 
the target prevents to use commom broad band filters. Individual frames were taken using the 
minimum possible exposure time of 1.2s. The total exposure time per filter was 30s in all cases. 
The data were reduced using the SOFI pipeline provided by ESO together with the Gasgano file 
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manager. Photometric calibration was obtained by comparison with field of view stars included in 
the 2MASS catalog. The zero points were derived using between 6 to 8 stars with accurate 2MASS 
magnitudes, the standard deviation was about 0.03 mag in all cases. 

In addition, SOFI was used to obtain a spectrum using the high resolution grism covering 
the K band with a slit width of 1". A total exposure time of 560 s was dedicated to the EX 
Lup spectrum, split into 8 individual frames of 70s each. A close solar type star (HIP78466) was 
observed to correct for telluric absorption. The data were also reduced using the SOFI pipeline. 

The SOFI and the GROND near- infrared photometries agree with each other within the un- 
certainties in the H and K bands. In the J band the difference is about 3.8 a, which indicate real 
difference, however with low significance. Such difference can be caused by short time-scale (several 
hours) low amplitude intrinsic photometric variability. 

The total spectral line contamina tion in EXors is usuall y weaker at near-infrared wavelengths 
compared to the optical (see also e.g. iLorenzetti et al.l 120071 ). To verify this in the case of EX Lup 
we used the SOFI K-band spectrum for the NB2.195 filter. Since for t he remaining two fi lters we 



did not have simultaneous spectroscopic data we used the spectra f rom (IKospal et al 



NB1.71 and NB1.215 filters, respectively. Although the data in (IKospal et al 



2011 



201 ih for the 



) were taken 

three months later than our photometry, the optical brightness of the source was about the same at 
both epochs and we assumed that the strength of the line emission was also similar on these dates. 
Since the observed line contamination was less than 1% in all three filters, which is significantly 
less than the photometric accuracy, we did not correct the photometries for line contamination. 



2.5. SPITZER IRS 



EX Lup was observed with the Infrared Spectrograph (IRS; iHouck et al.ll2004l ) onboard the 
Spitzer Space Telescope at four epochs; 18 March 2005 (PID: 3716), 21 April 2008 (PID: 477), 
10 October 2008 (PID: 524) and 7 April 2009 (PID: 524). The data of the first two epochs were 
published in [Abraham et al.l (j2009l ) and we refer to that paper for the details of the data reduction 
procedure of those spectra. 

At the last two epochs EX Lup was measured using only the low-resolution modules (R=60- 
120). Our spectra are based on the droopres products processed through the S15.3.0 version of the 
Spitzer data pipeline. First, the background has been subtracted using associated pairs of imaged 
spectra from the two nod positions, also eliminating stray light contamination and anomalous dark 
currents. Pixels flagged by the data pipeline as being "bad" were replaced with a value interpolated 
from an 8 pixel perimeter surrounding the flagged pixel. The spectra were extracted using a 6.0 
pixel and 5.0 pixel fixed-width aperture in the spatial dimension for the Short Low and the Long 
Low modules, res pectively. The low- level fringing at wavelengths >20;um was removed using the 
irsfringe package (jLahuis Sz BoogertI l2003l ) . The spectra were calibrated with a spectral response 
function derived from IRS spectra and MARCS stellar models for a suite of calibrators provided 
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by the Spitzer Science Centre. To remove any effect of pointing offsets, we matched orders based 
on the point spread function of the IRS instrument, correcting for possible flux losses. 



2.6. SPITZER MIPS 



We obtained photometry at 70 and a lo w resolution spect rum of EXLup, using the Multi- 
band Imaging Photometer for Spitzer (MIPS; iRieke et al.ll2004l ). The 70 //m imaging was per- 
formed in photometry mode. The data reduction was started with the Basic Calibrated Data 
(BCD) products generated by the pipeline (version S17.2) developed at the Spitzer Science Cen- 
ter. We perfo r med c olumn spatial filtering and time median filtering on the BCD images following 
Gordon et al.l (j2007l ). The filtered BCD data were coadded an d corrected for array distor tions with 
the SSC MOPEX (Mosaicking and Point source Extraction, iMakovoz fc Marleaul l2005l ) software. 
Outlier pixels were rejected using a 3a clipping threshold. The output mosaic was resampled to 
4"pixel~^. 

We used a modified version of the IDLPHOT routines to detect the source on the final image 
and to obtain aperture photometry with an aperture radius of 16" and with sky annul us between 18" 



and 3 9". The aperture correction, appropriate for a 60K blackbody, was taken from I Gordon et al 



(j2007l ). The final uncertainties were computed by quadratically adding the internal error and the 
absolute calibration uncertainty for which we adopted 7% (see MIPS Data Handbook). 

The low-resolution far-IR (55-95/im; R~ 15-25 ) spectrum of EXLup was obtained in the 
spectral energy distribution (SED) mode of MIPS. Six observing cycles were performed, each con- 
sisting of 10 s long on- and off-source exposures. The on-source and the background positions were 
separated by 1". We began our data reduction with the BCD images (pipeline version S17.2) and 
MOPEX was used to perform further processing steps (combination of data, background removal, 
application of the dispersion solution). The spe ctrum was extr acted in an 5 pixel wide aperture. 
The aperture correction factors were taken from lLu et al.l (|2008l ) . 



2.7. VLTI MIDI 



EXLup w as observed with MIDI (Leinert et al.ll2003l ). the N-band instrument on the VLTI 



interferometer (|Glindemann et al.M2000l ). Spectrally resolved interferometric visibilities were ob- 
tained between 8 and 13 fim with a spectral resolution of i2~35. Observations were performed 
at two epochs, during the nights starting 21 June 2008 (baseline U2U3, projected length 37.4 m, 
position angle 57.6° E of N) and 16 July 2008 (baseline U1U4, projected length 121.1m, position 
angle 73.4° E of N). Data were obtained in "high sense" mode, i.e. each interferometric observation 
in which the beams coming from both telescopes are combined coherently is followed by photomet- 
ric observations in which first only the light coming from telescope A, then from telescope B, is 
measured, and the intensities of the individual beams are determined. 
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Calibration stars were observed to determine the system response and atmospheric trans- 
parency. During the night of 21 June 2008 the cahbrator HD 149447 (K6III, 6'=4.68±0.05 mas, 



Cohen et al.lll999l ) was observed immediately following the EXLup observation. During the night 



of 16 July 2008 the same calibrator was observed immediately before and after the EX Lup obser- 
vations. We used the EWS package to reduce the data, and verified that reductions using MIA 
yield consistent results. 



2.8. VLT VISIR 

Low-resolution (R~250) N-band spectroscopy of EXLup was performed with VISIR at the 
VLT. We observed the source during two epochs, first around 25 July 2008 (for which data from 
the nights of the 24, 25, and 27 of July were combined), and second in the night of 28 August 
2008. Standard chopping and nodding techniques were applied to correct for the high celestial and 
instrumental background inherent to mid-IR observations. 

During the July observations, HD 149447 was observed as telluric and flux calibrator immedi- 
ately following each observation of EXLup, at essentially identical airmass. Du ring the observation s 



in August we used HD 196171. Both calibrators are taken from the catalogue of I Cohen et al.l (| 19991 ). 
that provides spectral types and angular diameters for each star, which we used to compute Kurucz 
model spectra to calibrate the absolute flux levels. 



2.9. APEX LABOCA 

EXLup was observed at 870/im using the 295-element LABOCA (Large Apex Bolometer Cam- 



era, Siringo ( 


it al. 


Glisten et al. 


2006) 



tion mode using a spiral pattern. This observing mode provides a fully sampled map of the total 
field-of-view of LABOCA. The total on-source observing time was 2.2 h. The data processing was 
performed by the Bolometer array Analysis Software (BoA) package following the data reduction 
steps (correction for the atmospheric opacity, flux calibration, removal of correlated noise using 
an iterative method, removal of spikes) outlined in the BoA User Manual (version 3.1). For the 
flux calibration we used the sources IRAS16293, B13134, G34.3.The co-added final map had pixels 
with sizes of 6". Aperture photometry was used to extract flux from the map using the position of 
EXLup. 
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3. Results 
3.1. Light curve 



The brightest ever recorded outburst of EX Lup was announced in February 2008 (|Jonesll2008l ) . 
The hght curve, composed from visual estimations and V-band observations, is presented in Fig.[TJ 
Unfortunately the beginning of the outburst was not observed, one can only extrapolate for the 
starting date. The last optical measurement (ASAS) in the quiescent phase was on 10th Oct 2007 
(~110 days before the peak brightness), when brightness was V=12.7 mag. We estimated the total 
length of the outburst to be about 300 days, which is an average value within the class. 

The brightness of EX Lup at the peak of the o utburst was about V=8 mag. Taking the quiescent 



brightness to be V~ 13 mag (jSipos et al.l |2009| ) one gets a V-band amplitude of 5 mag for the 



outburst. With this value, the 2008 outburst of EX Lup belongs to the strongest ones observed in 
EXors. This is an unusually high value compared also to the historical record of EX Lup. The only 
comparable outburst was observed in the 1950s. 

The shape of the light curve during the 2008 outburst of EX Lup showed some peculiar features. 
The observations began about three weeks before the peak of the outburst. The brightness of the 
source rose steeply until it achieved a peak brightness of V=8 mag. The star remained this bright 
only for a few days, then it faded to about 9.5 mag. Then the average brightness decreased by 
another 2 mag in the next 200 days. One of the most interesting properties of the light curve was 
the quasi-periodic brightness variation seen on top of the plateau-like slow fading. The amplitude 
of these fluctuations increased with time, while their period remained roughly constant, about 35 
days. 



3.2. Spectral Energy Distribution 

The SEDs of EX Lup on 21 April 2008 and in the quiescent phase are compared in Fig.[2j 
The outburst-to-quiescent flux ratio was the highest in the U-band and it decreased gradually 
with wavelength. Although EX Lup was already fainter on 21 Apr 2008 than at the peak of the 
outburst, at optical wavelengths it was still more than an order of magnitude brighter in the outburst 
compared to the quiescent phase. The shape of the SED in the optical domain is consistent with 
blackbody emission at ~ 6500 K, which is far higher t han the effective temperature of the star in 



quiescent phase (3800 K. I Gr as- Velazquez &: Rayll2005l ). Although the absolute flux levels in the 



GROND i and z filters are somewhat higher than what would be expected for the aforementioned 
blackbody emission, the i-z color fits well to the modeled emission. The discrepancy in the absolute 
flux level can be explained by the relatively high calibration uncertainty in the GROND optical 
bands. In the near- infrared domain the flux increased only by a factor of about 2.5 in the outburst 
compared to the quiescent phase. Moving to longer wavelengths the ratio of the outburst-to- 
quiescent fluxes increased again to about an order of magnitude at 5.5 ^m, which is the shortest 
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wavelength of the Spitzer IRS spectra. Longwards of 5.5 fim the ratio of outburst-to-quiescent 
fluxes decreased with wavelength. 

The brightening was certainly caused by the increased accretion rate onto the star during 
the outburst, indicated by the increasing luminosity in the accretion indicator lines. In Fig. [3] 
we c ompare the s pectra l region around the Br 7 line in the outburst and in the quiescent phase 



from ISipos et al.l (j2009l ). It can be seen that the luminosity of the Br 7 was far h igher in the 



outburst than in the quiescent phase. Using the relationship of lMuzerolle et al.l (|l998l ) we derived 



an accretion rate of 2 • 10~^MQ/yr in the outburslj^te, that the hydrogen recombination lines are 
not only indicators for accretion, but also mass loss. Therefore, the derived accretion rates should 
be considered as upper limit. This is ab o ut th ree orders of magnitude higher than the value in 



the quiescent phase found by ISipos et al.l (|2009l ). We want to note, that different lines were used 
to derive the accretion rate in the outburst (Br 7) and in the quiescent (Pa/3), since none of the 
accretion indicator lines were measured both in the quiescent and in the outburst. The usage of 
different accretion indicators might introduce some additional uncertainty when a comparison of 
the accretion in the outburst and in the quiescent is done. This uncertainty, however, does not 
affect our conclusion that the 2008 brightening of EX Lup was certainly caused by the increased 
accretion rate onto the star during the outburst, indicated by the increasing luminosity in the 
accretion indicator lines. 



3.3. Mid-infrared spectra 

In this section we will describe the mid-infrared spectra measured towards EX Lup in our 
observing campaign. In order to avoid confusion, we will use the following naming convention. The 
ensemble of features between 8 and 13 ^m will be called the " 10 ^m feature complex" . The " 10 fim 
feature" term will be used for the forsterite band peaking at 9.9 ^m. The Spitzer IRS spectra of 
EX Lup are presented in Fig.Hl In order to show also the weaker dust features, a spline continuum 
was fitted to and then subtracted from each spectrum. The continuum subtracted spectra are 
shown in Fig.[5j 



The pre-outburst spectrum from March 2005 was already analyzed by ISipos et al.l (|2009l ). 
They showed that the dust in the disk atmosphere of EX Lup is dominated by small amorphous 
grains with a mass- weighted average grain size of 0.5 //m. They also showed that the abundance of 
crystalline silicates in the disk atmosphere is negligible. 

The Spitzer IRS spectrum from the outburst, taken on 21 April 2008, showed remarkable 
changes. First of all, the absolute flux level increased by about an order of magnitude at the 
shortest wavelength of the spectra (5.5 fim) and by about a factor of five at the longest wavelength 
(~35^m). The amount of brightening decreased with wavelength shortwards of about 20 ^m, but 



2n 
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i t was completely wavelength independent longwards of 20 /xm. As reported by lAbraham et al 



(|2009l ) the spectrum showed a small, but narrow and sharp feature at 10 /xm on top of the 10 fim 



feature complex (see Fig.[5]). A shoulder at 11.2 //m and another feature at 16 u m also appeared 



i n the outburst spectrum which were not present in the pre-outburst spectrum. [Abraham et al 



(|2009l ) identified the source of these new bands as crystalline forsterite. It was also reported that 



no crystalline feature was seen in the outburst spectrum longwards of 20 /im. 

The third Spitzer IRS spectrum, taken on 10 October after the end of the outburst, showed 
further changes. As can be seen in Fig. HI the absolute flux level decreased compared to the outburst 
spectrum. The flux level was still higher than in the quiescent spectrum by about a factor of 1.7, but 
this factor was wavelength independent in the whole Spitzer IRS wavelength domain. The shape 
of the spectral features between April 2008 and October 2008 also changed. In Fig. [5] it can be seen 
that the 10 micron feature complex became broader and the 11.2 ^um shoulder became stronger, 
relative to the underlying broad amorphous feature, in October compared to the April spectrum. 
The most remarkable changes occured, however, at longer wavelengths. Apart from the narrow 
16 fim band of forsterite, new bands appeared in the October 2008 spectrum at approximately 
23 //m, 28 /im and 33 /xm. These wavelengths corresp ond to the peak w avelengths of the strongest 



forsterite bands in this wavelength interval (see e.g. iKoike et al.l 120031 for a list of forsterite band 



positions). Thus we identified the source of these bands to be crystalline forsterite. 

The last Spitzer IRS spectrum was taken on 7 April 2009. The absolute flux level decreased 
further compared to the first post-outburst spectrum from about five months earlier. Longwards of 
20 /xm the flux decrease was independent of the the wavelength. The flux variation showed, however, 
clear wavelength dependence shortwards of approximately 20 /xm, where the brightness decreased 
the most. The flux at 5.5 /xm decreased to about half of the pre-outburst value, while longwards of 
20 /xm the flux in the last IRS spectrum was approximately 30 % above the pre-outburst value. The 
shape of the spectral features longwards of 13 /xm remained unchanged within the uncertainties (see 
Fig. [5]), while the 10 /xm silicate feature complex clearly changed compared to the previous epoch. 
It can be seen in Fig.[5l that the 10 /im forsterite feature on top of the broad amorphous feature was 
still present, but the shoulder at 11.2 /xm became much weaker in April 2009 than it was in October 
2008. The general shape of the 10 /xm feature seems to be somewhere between the outburst and 
pre-outburst spectra. 

Due to the difference in the noise levels between the outburst and the post-outburst spectra, 
we cannot exclude the possibility from the Spitzer IRS spectra alone that weak crystalline features 
(of the same strength as in the post outburst spectra) longwards of 20um were already present in 
the outburst spectrum. In the post-outburst spectra the peak of the 24 and 28 um features is 9 
and 6 times the corresponding flux uncertainty above the underlying continuum, respectively. The 
absolute noise level of the outburst spectrum is, however, worse than the post outburst spectra by 
a factor of about 3.4. If we measured the post-outburst spectra with the same noise level as we 
had in the outburst spectrum, the detection of the 24 and 28 um features would be about 2.6 and 
1.8 times the corresponding flux uncertainty, respectively, i.e. they could be hidden by the noise. 
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We will further discuss this question in Sec. l4.2.r] 

In the Spitzer IRS spectra, the disk of EXLup remains unresolved allowing us to investigate 
the integrated spectrum of the disk as a whole only. To obtain information about the spatial 
distribution of the newly formed crystals we used MIDI on the VLTI. In our first measurement in 
the U2U3 configuration (22 Jun 2008) with a projected baseline length of 37.4 m EX Lup remained 
spatially unresolved. The second measurement in the U1U4 configuration on 16 July 2008 (projected 
baseline lengt h 121.2 m), howe v er, re solved the source. In the following we will use the naming 



convention of Ivan Boekel et al.l (j2004l l and we will refer to the correlated spectrum as "inner disk 
spectrum" . We also subtracted the correlated spectrum from the total spectrum and the resulting 
uncorrelated emission will be called the "outer disk spectrum". In Fig.^Left we present the inner 
disk spectrum of EXLup. The spectrum peaks at 11.3 fj,m which corresponds to the wavelength of 
the strongest forsterite peak in the 8-13 fim region. The inner disk spectrum of EXLup resembles 



that of HD 142527 which is known to originate in regions with > 90 % crystallinity (Ivan Boekel et al 



20041 ). suggesting also very high crystallinity in the EXLup inner disk spectrum. The outer disk 
spectrum of EXLup is shown in Fig.MRight. It is broader compared to the pre-outburst total 
spectrum and peaks at about 10.2 ^m. The feature is smooth and lacks any substructure, which 
would be the signature of crystalline silicates. Interesti ngly, the outer disk spectrum of EXLup 



is very similar to the 10 //m feature complex of FUOri (jQuanz et al.l 120071 ). iQuanz et al.l (|2007l ) 
have shown that the 10 /xm feature of FU Ori is known to emerge from a region where significant 
grain growth has taken place, indicating also the presence of large (few /um) grains in the outer 
disk atmosphere of EXLup. 

Two spatially unresolved spectra were taken with the VISIR instrument on VLT on 25 Jul 
2008 and on 28 Aug 2008 (see Fig.[TTJ. Although the VISIR spectra have better signal-to-noise 
ratio than those taken by the MIDI instrument, the quality of the spectra is not as good as that of 
Spitzer IRS, allowing only moderate analysis of the spectral features. In July 2008, at the epoch of 
the first VISIR measurement the source was at about the same brightness level as in Apr 2008, when 
the outburst Spitzer IRS spectra was taken. The absolute flux level of the two spectra agrees well 
with each other, but the peak-to-continuum ratio is lower in the VISIR spectrum than than in the 
Spitzer IRS (1.6 vs. 1.4). The 10 fim feature complex seems to be broader in the VISIR spectrum, 
too. The broader and shallower features are usually interpreted as indication for the presence of 
larger, micron-sized grains. Similar to the MIDI data the VISIR spectrum points towards the 
presence of large grains, too. The Aug 2008 VISIR spectrum was taken at about the same optical 
brightness level as the first Spitzer IRS post-spectrum. Due to the lower signal-to-noise ratio of the 
VISIR spectrum no significant diff'erence can be found in the shape or the strength of the Spitzer 
and VISIR spectra above the noise level. 
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Discussion 



4.1. Crystal formation 



Abraham et alj ( 2009 ) reported the in situ formation of forsterite crystals via annealing in the 
disk atmosphere of EXLup during its 2008 outburst. The crystalline features we found in the 20- 
35 jum wavelength i nterv al are also associated with crystalline forsterite, supporting the conclusion 
of Abraham et al. ( 2009 ). It is interesting that no features of other crystalline dust species, apart 
from forsterite, can be seen in the mid-infrared spectra. For instance, enstatite (the magnesium- 
end member o f the pyroxene solid so l ution series), which is frequently obser ved in protoplanetary 
disks (see e.g.. Ivan Boekel et al.ll2005l . I Watson et al.ll2009l . lJuhasz et al.ll20ld ). is not present in the 
spectra of EXLup. Laboratory annealing experiments of amorphous silicates show that forsterite 
will form independe ntly of the starting stoichiometry if the grains are small (< 1 /um) and porous 
( Fabian et al.ll200d ). For larger (few ^um) and compact particles, annealing of amorphous grains 
with pyroxene stoichiom etry results in the formation of enstatite. Using the activation energies 
from ([Fabian et al.ll2000l ) the annealing time of enstatite in the 1000-1300 K temperature interval is 
about an order of magnitude longer than that of forsterite, but it decreases rapidly with increasing 
temperature. At 1200 K the calculated annealing time of enstatite is on the order of a minute. 
It is unlikely that dust grains in the outburst are heated up to 1200 K then cooled down below 
1000 K on a one-minute time-scale to prevent the formation of enstatite. The only possibility for 
the formation of forsterite without enstatite production is that the parent amorphous grains were 
sub-micron sized porous grains. Indeed, the analysis of the Spitzer IRS spectrum of EXLup from 
the qui escent phase show the presence of amorhphous grains with a mass-averaged grain size of 
0.5 fim (|SiDos et ahlbood ). 



All crystalline silicate features, so far identified in protoplanetary disks, were associated with 
iron- po or, magnesium rich crysta ls (see e.g. van Boekel et al1l2005 : Bouwman et al ] l2008l :l ISargent et al 
20091 or lHenning &: Meeua 120061 for a review). An important consequence of the absence of iron is 
the dramatic decrease of the optical and near-infrared opacity of the grains. Iron-poor silicates are, 
therefore, cooler than ferromagnesian silicates at the same distance from the central star. Compar- 
ing the positions of the obs erved forterite featu res in the spectrum of EX Lup with those of synthetic 
and natural olivines from iKoike et al.l (|2003l ). we estimated the iron content of the lattice to be 
less than 10 % with respect to magnesium. Such iron-poor crystals can be formed from iron-poor 
amorphous silicates, but also from ferromagnesian silicates if the annealing occurred under reducing 
conditions. Laboratory experiments showed that in the latter case the redu ced iron will remain i n 
thermal contact with the resulting crystals in the form of metallic inclusions (iDavoisne et al.ll2006l ). 
Such iron inclusions do not affect the positions of the forsterite bands while the grain as a whole 
will have higher optical to near-infrared opacity than that of isolated, pure forsterite crystals (see 
e.g. 



Ossenkopf et al.lll992l ). 



Besides iron, amorphous carbon is the other important opacity source in protoplanetary disks. 
Carbon, as well as iron, has a featureless opacity curve and high optical to near-infrared opacity 
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(j Jaeger et al.lll998l ). ISipos et al.l (|2009l ) used amorphous carbon (20% in terms of the total dust 
mass) in their radiative transfer modehng of the quiescent phase SED of EX Lup. Amo rphous 



carbon, however, is not stable above 1000 K as it oxidizes quickly and forms CO (|Gailll200ll ). Since 



annealing is expected to occur above 1000 K, we do not expect amorphous carbon to be in thermal 
contact with silicate crystals. 

Based on the above constraints, we propose two scenarios to be tested for the formation of 
forsterite crystals in the disk of EXLup. In scenario A the parent grains were porous ferromagnesian 
silicates in thermal contact with carbon. As the outburst set in and the temperature rose above 
1000 K, the silicate grains crystallized. The formation of CO from amorphous carbon lowered 
the partial pressure of oxygen providing reducing conditions under which iron-poor crystals can 
form from ferromagnesian silicates. To explain the observed variations of the strength of forsterite 
features an additional mechanism, e.g. radial/vertical mixing is required. 

In scenario B the parent amorphous silicate grains were iron-poor sub-micron sized porous 
grains in thermal contact with carbon. As the temperature rose above 1000 K at the beginning of 
the outburst, the silicate grains were annealed to crystalline forsterite, while the amorphous carbon 
oxidized and formed CO. Due to the low opacity of the crystals (in the absence of carbon) their 
temperature decreased. This resulted in the weakening of the forsterite features in the 10 fim region 
while in the same time the features in the 20-35 /im region became stronger, which might explain 
the observed spectral variations. 



4.2. Modelling 



We modeled the SED, the mid-infrared spectra and the interferometric visibilities sin iultane- 
ously using the 2D radiative t ransfer (RT) code RADMC (Dullemond DominikI l2004al ) and a 
1-l-lD turbulent mixing code (jDullemond &: DominikI l2004bl ) . The goal of our modeling was to 
find answers to the following three questions. Are the freshly produced crystals isolated iron-free 
particles or are they in thermal contact with iron? Which physical process was responsible for the 
observed variation of the mid-infrared solid-state features? What kind of constraints can we put 
on the mechanism of the outburst? We calculated, therefore, seven series of models to simulate the 
two crystal formation scenarios (series A-, and B-), with (-AD-) and without (-NAD-) an optically 
thick accretion disk and with different strength of the vertical turbulent mixing (-aO.Ol, -aO.l, 
-al.O). Details of the modeling procedure are described in Appendix A. We now briefly describe 
the modeling strategy in the different model series. 



Series A-AD-al.O, A-AD-aO.l, A-AD-aO.Ol. First, we started with a disk containing 
only amorphous silicates an d amo rphous carbon. The parameters of the starting disk model 
were taken from lSipos et al.l ((20091) and are sumn iarized in Tab.Bl The inner radius of the disk 
was taken to be 0.3 AU (instead of 0.2 AU as in lSipos et al.ll2009l ) to get a better agreement 
between the model and the quiescent phase Spitzer IRS spectrum in the 5-10 /xm wavelength 
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interval. After the accretion rate has been set, we ran the RT code to obtain the temperature 
structure at the beginning of the outburst. Amorphous dust grains were replaced by forsterite 
grains with iron inclusions wherever the temperature in the disk exceeded 1000 K. Then the 
turbulent mixing code was run to simulate vertical mixing over a 10-day period. Turbulent 
mixing was treated as a diffusion process using the a prescription ([Shakura &: Sunvaevlll973l . 
see also Appendix A) for the diffusion coefficient. The strength of the turbulence was con- 
trolled via the value of a, for which we used 0.01, 0.1 and 1.0. The last step was to update 
the density structure according to the results of the mixing code. The whole modeling cycle 
was then repeated. Accretion heating was simulated with an optically thick accretion disk as 
well hot spot on the stellar surface. 

Series A-NAD-al.O, A-NAD-aO.l, A-NAD-aO.Ol. Same as the A-AD- models, but no 
accretion disk was present in the models. The total accretion luminosity was assumed to be 
radiated away from the hot spot. 

Series B-AD. We started the modeling in the same way as in the A-AD- models. Amorphous 
dust grains were replaced with isolated iron-free forsterite crystals wherever the temperature 
in the disk exceeded 1000 K. Amorphous carbon was also removed from the dust phase above 
1000 K. In this model series we did not use the vertical mixing code, we only updated the 
accretion rate and calculated models on each 10th day. 



For the amorphous grains we applied the dust model of ISipos et al.l (|2009l ). which was used 
to model the quiescent phase SED. During the modeling we used two grain models for forsterite 
crystals. For isolated iron-free forsterite crystals we used mass absor ption coefficients calculated 
from the optical constants using Distrib ution of Hollow Sph eres (DHS. lMin et al.ll2003l ) scattering 
theory. We used t he optical constants o f ISogawa et al.l (j200g) longwards of 2 //m and supplemented 
them with that of lScott &: DulevI (|l996l ) shortwards of 2 /im. For composite grains, where forsterite 
and iron are in thermal contact with each other, we calculated the optical constants of the aggregate 



froni those of the constituents first, using effective medium theory (APMR mixing rule, iMin et al 



20081 ) . For the constituents of the aggregate we used the DHS polarizability. Finally Mie theory was 
applied to calculate the mass absorption coefficient of the aggregate. The same optical constants 
w ere used for f orster ite and for the isolated crystals, while for iron we used the optical constants 
of lOrdal et al.l (119881 ) between 0.66 /im and 285 ;um and extended it with the optical constants of 



Pollack et al 



(|l994l ) shortwards of 0.66 ;um and longwards of 285 /xm. 



In Fig.[7| we compared the observed SED of EXLup in the outburst to model calculations. 
Although our models reproduced the general shape of the observed SED, there were quantitative 
differences. The largest differences between our models and the observations occured in the 1- 
8 [xm. wavelength interval (A-AD models) and in longwards of 60 [iva.. The former is related to the 
presence or absence of an optically thick accretion disk inwards of 0.3 AU. We will discuss this 
problem in Sec. 14. 2. 21 
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The most likely explanation why our models overestimated the observed fluxes longwards of 
60 fim is that the disk of EX Lup was not in an equilibrium state. RADMC, as well as other RT codes 
used to model protoplanetary disks, calculates the temperature structure of the disk in radiative 
equilibrium. The emission longwards of 60 /im is dominated by the disk interior which is optically 
thick to the heating radiation of the star and the disk atmosphere above it. The high optical 
thickness prevents the disk interior to react to the variable irradiation luminosity quickly. Since 
the optical depth in the far-infrared to millimeter wavelengths is lower than at shorter wavelengths 
the observer can see deeper into the disk observing those regions, that are still not accommodated 
to the changed irradiation luminosity. The optically thin, upper layers of the disk can, however, 
react to the variable irradiation luminosity almost immediately. 



4-2.1. Radial and vertical mixing of solids 

The first question we tried to answer with our modeling was which of the two crystal forma- 
tion scenarios, outlined in Sec. 14. 11 was at work in the disk of EX Lup. In Fig. [8] we compared the 
calculated mid-infrared spectra of models A-Aclfl and B-AD. As can be seen in Fig.[8ti-c, the posi- 
tions of forsterite bands in the observed and modelled spectra matched very well in models A-AD. 
The strengths of the forsterite peaks are dependent, however, on the strength of the vertical mix- 
ing. The positions of the forsterite peaks in the B-AD models were shifted systematically towards 
shorter wavelengths compared to the observations (see Fig.[8]i). The shift of the peak position was 
certainly related to temperature/radiative transfer effects since in the mass absorption coefficients 
there was no difference in the peak positions of isolated forsterite crystals and forsterite-iron ag- 
gregates. Moreover spectra of B-AD models showed also stronger forsterite features, especially in 
the post outburst epochs, compared to the observed spectra. Thus, our modeling suggests that the 
freshly produced forsterite crystals had to be in thermal contact with iron, which excludes scenario 
B as a crystal formation mechanism. 

The next question we addressed was whether vertical mixing can explain the observed variation 
in the shape of the mid-infrared solid-state emission features. Without considering vertical mixing, 
the apparent crystallinity in spatially unresolved spectra depends on the ratio between the size of 
the crystallized region and the total (crystalline -|- amorphous) silicate emission zone seen at 10 /im. 
The former depends on the irradiation luminosity at the peak of the outburst and does not change 
with time. The total silicate emission zone, ho wever, scales linearly with the actual irradiation 



luminosity (see, e.g., iKessler-Silacci et al.l 120071 ) . If the irradiation decreases, the total silicate 



emission would shrink accordingly. Since the size of the crystalline region remains unchanged, 
the apparent crystallinity increases. So far we assumed that in the crystalline region the value of 
crystallinity is 1.0, i.e., every dust grain is in crystalline form and this does not change with time. 



^The spectra of models A-AD and A-NAD are very similar to each other, therefore, to avoid confusion we show 
only the spectra of models A-AD. 
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The effect of vertical mixing is to decrease the crystallinity with time in the crystalhne zone by 
stirring up amorphous grains from deeper layers of the disk and, at the same time, mixing down 
crystals. In this way vertical mixing might compensate the effect of the shrinking total silicate 
emission zone. 

In Fig.[8ti.-c we compared the observed spectra in the 7.5-13.5 fim region with those of models A- 
AD with different values of a. In our weakest mixing case (a = 0.01, Fig. [8^), our models predicted 
stronger forsterite features than observed. In contrast, models with the strongest possible mixing 
(a = 1.0, Fig. [8^), predicted too weak crystalline features compared to our observations. The best 
agreement between our models and the observed spectral features in the 7.5-13.5 ^m was achieved 
using a = 0.1. The strength of the 10 fj,m feature complex is well reproduced for the outburst (Apr 
2008) and the second post-outburst spectra (Apr 2009), but our model overestimates the feature 
strength compared to the spectrum taken in Oct 2008. The disagreement between the model and 
the observed spectra is about 35%. Since the emission at the 10 fim is optically thin it reacts to the 
changing irradiation immediately and the infrared radiation is proportional to the irradiation flux. 
This means that the 35% difference between our models and the observed flux level translates to 
a difference of ~0.3 mag difference in the optical. Since we had no optical photometry at the time 
when the Spitzer IRS spectrum was taken, only a few days after, the difference could be caused 
by the extrapolation to the optical brightness used in our modeling. In Fig[9^-c we compared our 
models to the observed spectra in the 13.5-36 ^m wavelength interval. Although our models with 
a = 0.1 matched the observed spectra in the 7.5-13.5 /im region well, the observed forsterite features 
longwards of 20 /xm are far stronger than in any of our models independently of the value of a (see 
Fig.ED. 

Two explanations could be used to solve the discrepancy between the observed and modeled 
crystalline features longwards of 20 /im. A solution could be that, in reality, the size of the crys- 
tallized region was higher than in our models. This can be due to a higher peak luminosity of 
EXLup in the outburst or due to a lower crystallization temperature. In this case, the outburst 
spectrum could have contained some weak forsterite emission features longwards of 20 /im, which 
were hidden by the noise in the measured spectra. The uncorrelated MIDI spectrum, however, 
excludes this possibility. We tested this scenario with our 2D RT model. If the crystallized region 
is larger than about 1.5-1.8 AU, the calculated uncorrelated spectra of our models show strong 
forsterite bands, which are not observed in the MIDI spectra. Moreover, models which reproduced 
the strength of the forsterite features longwards of 20 /im in the spatially unresolved post-outburst 
spectra, overestimated the strength of the crystalline features in the 10 //m region. Using our mod- 
eling framework with vertical mixing, we were not able to reproduce the strength of the crystalline 
features both shortwards and longwards of 20 /im in the post-outburst spectra, and the spatially 
resolved mid-infrared interferometry simultaneously. 

The other explanation would be that forsterite crystals were not present in the disk atmosphere 
outwards of 1.5-1.8 AU in July 2008, but they were transported to a few AU by October 2008. The 
transport of crystals over a 2 AU distance within three months implies an average radial velocity of 
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about 38km/s. Radial mixing is far too slow for such a transport velocity. Instead, it is likely that 
in this scenario the crystals are transported by e.g. a wind driven transport mechanism. Indeed, 
optical and near-infr ared spectroscopic observations suggest the presence of a wind during the 
outburst of EXLup ( lAspin et al.ll2010l : iGoto et al.ll201ll : iKospal et al.ll201ll . Sicilia-Aguilar et al., 
in prep). The quantitative modeling of this scenario is, however, beyond the scope of this paper. 



^.2.2. Outburst mechanism 
EXor outbursts are usually considered to be a 'down-scaled' version of FUor outbursts, sug- 



gesti ng the same triggering mechanism of the outbursts in the two classes (jHartmann &: Kenvon 



1996). The most wid ely accepted model of FUor outbursts is the self-regulated thermal instability 
(TI. iBell Linlll994l ). Although being the m ost widely accep ted model, even TI cannot explain 



all observational properties of FUors (see e.g. IZhu et al.l 120071 ). In case of TI the released viscous 



energy in the disk midplane is trapped in the disk above a certain threshold surface density value. 
The radiative trapping is caused by the opacity of H~ ions in the partially ionized disk regions. The 
effective viscosity as well as the mass transfer rate in the disk within the ionized radius increases. 
The ionization front propagates outwards in the disk up to a certain radius then retreats again. 
The gas inwards of the ionization front is accreted onto the central star. 

In the following we will investigate the TI model in the context of our observations. TI 
requires an optically thick accretion disk (due to the presence of H" ions) whi ch domina.tes th e 



SED of FU Ori-type objects at optical and near infrared wavelengths (see e.g. IZhu et al.l 120071 ). 
The comparison of our A- AD and A-NAD models to the observed SED of EXLup showed that 
models without an optically thick accretion disk fit the SED better in the 1-8 //m wavelength range 
(Fig. [7]). In order to quantify the goodness of the SED fit we calculated and compared the of the 
fit for all model series. The number of tuned parameters were only two, the value of a (except for 
B-AD models) and whether or not an optically thick accretion disk was present within the inner 
hole of the dust disk.. The number of fitted data points were 1435. The values for the B-AD 
model was 3409. The A-AD-aO.Ol, A-AD-aO.l and A-AD-al.O models resulted in a of 1423, 
1372,1290, respectively. Finally the A-NAD-aO.Ol, A-NAD-aO.l and A-NAD-al.O models provided 
a x^ of 533, 470 and 492, respectively. It can be easily seen that the goodness of fit is very sensitive 
to the presence (or rather to the lack) of an optically thick accretion disk, while only marginal 
dependence is visible on the value of a. Independently of the value of the viscosity parameter, the 
best fit model family is the A-NAD, while the absolute best fit model is A-NAD-aO.l Mid-infrared 
visibilities of models without an optically thick accretion disk gave also better agreement with the 
observations (see Fig. llOp . This suggests that there was either no optically thick accretion disk 
present during the outburst of EXLup, or the emission at optical wavelengths came from a narrow 
ring of hot gas in the inner edge of the gas disk. 



One of the main differences between FUor and EXor outbursts is the time-scale over which 
the outburst occurs. The TI model alone, without any additional perturbation, is not able to 
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reproduce the rise-time in the case of FUOri and of V1515 Cygni (jClarke et al.ll2005l ). In these 
cases the brightening happened over about a year, which is an order of magnitude longer than in 
the case of EXLup (few weeks). The relatively short duration time of the outburst of EXLup 
puts constraint also on the outer radius, from where gas can be accreted onto the star. The mass 
transfer in the disk occurs by viscous diffusion, the time-scale of which can be estimated using 



7 

D 



an' 



(1) 



where D is the diffusion coefficient {D = aCsHp), r is the radius, Hp is the pressure scale height 
of the gas and is the Kepler frequency. We used Eq.[T] to estimate the outer radius from which 
material can be accreted within 10 months. We assumed that r / Hp = 10 and for a we used the 
highest theoretically possible value, a = 1.0. The resultin g radius is ~ 0.1 2 AU, which is well 
within the fitted inner radius of the dust disk (0.2-0.35 AU, ISipos et al.ll2009l ). This implies that 
the reservoir of material for the 2008 outburst should have been present in the inner hole of the 
dust disk. This calculation implies that the hole in the disk, observed in the quiescent phase SED, 
should not be empty, but it is filled up with gas, which is optically thin in the continuum. On the 
other hand, if thermal instability takes place in the disk, the disk should be optically thick due to 
the presence of H~ ions, which should be visible in the SED. 

The mass of gas accreted onto the star during this outburst was about 10"'' Mq, and the 
calculation above implies that all of this gas should be located within 0.12 AU from the central 
star. If we assume that such high amplitude outburst, like that in the 1950s or in 2008, consume 
all the available gas within 0.12 AU a mass transfer rate of 2 • 10~^ Mq/ji or higher is required to 
maintain the observed amplitude and frequency (one in every 50 years) of these large eruptions. The 
estimated mass transfer rate is a lower limit, as we neglected all the smaller amplitude outburst, 
which occu r red s ince the 1950s, but it is more than an order of magnitude higher than what 
Sipos et al. (j2009l ) estimated. The req uired mass transf er rate in the disk for TI to operate was 



found to be 5 • 10~''MQ/yr or higher ( Bell Sz Linlll994l ). which is far higher than the estimated 
accretion rate in the quiescent phase of EXLup. If the outbursting region is fed with a lower 
rate, thermal instability does n ot occur, th e accr etion onto the star is stable and steady. The 
equilibrium curves calculated bv lBell &: LinI (jl994l ) also suggests that between 5 • lO~''M0/yr and 
about lO^^MQ/yr there is no stable solution for the accretion. If thermal instability sets in, the 
disk quickly flips to the outburst phase and the accretion rate increases to above lO~^M0/yr. 

The mass transfer rate in the disk is, however, not necessarily the same as the accretion rate 
onto the star. The accretion rate onto the star can be lower than the mass transfer rate in the 
disk due to, e.g., interactio n between the stellar magnetic field and the accreting gas (see, e.g.. 



Romanova et al.ll2002l . |2004| . D'Angelo et al. in prep), which causes the gas to pile up in the disk. 



The gas lines, used as accretion indicator, a re formed in the funne l flow between the inner radius of 
the gas disk and the surface of the star (e.g. iMuzerolle et al.lll998l ). therefore they only measure the 
accretion rate onto the star. Thus, the mass transfer rate in the disk can be significant even if the 
accretion rate derived from, e.g., the Br 7 line is zero. In FUor outbursts, all gas piled up over 10^ 
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years (estimated time between two subsequent FUor outbursts) is accreted onto the star since the 
accretion onto the star is as high as lO~^M0/yr. In contrast, the accretion rate of EXLup in the 
outburst was three orders of magnitude lower. If the mass transfer rate in the disk of EX Lup was 
5 • 10~'''Mo/yr (required for the TI to operate), EXor outbursts cannot accrete ah piled up mass 
from the disk onto the star. In this case gas accumulates continuously in the inner regions of the 
disk and at some points the system should go into a FUor-type outburst, meaning that EXors are in 
fact FUors in the "quiescent" phase. Since EXLup has a disk mass of 0.02 Mq and no evidence for 
a massive envelope has been found so far, the disk of EXLup would be entirely consumed by only 
a very few FUor eruptions (assuming M=10~^MQ/yr and an outburst duration of 50 yr), which 
makes this scenario unlikely. Our study suggests that the 2008 outburst of EXLup is caused by 
increased accretion onto the central star as indicated by the increase of the Br 7 line luminosity in 
the outburst compared to the quiescent phase. On the basis of the observational data we conclude, 
however, that it is unlikely that the self-regulated thermal instability is the responsible mechanism 
triggering the outburst. 



The recent work of lD'Angelo fc SpruitI (I2OIOI ) could provide an alternative explanation for the 
outburst mechanism of EXLup. If the magnetic truncation radius of a strongly magnetized star 
is close to the co-rotation radius the accretion onto the star will not be continuous anymore, but 
becomes episodic. After the end of an outburst, in the quiescent phase, the inner radius of the disk 
is located just outside of the co-rotation radius. The accretion onto the star is low and material 
piles up at the inner edge of the disk. The inner radius of the disk moves inwards and when 
it crosses the co-rotation radius a burst of accretion occurs that empties the piled up reservoir 
of material and moves the inner edge of the disk outside of the co-rotation radius again. The 
interesting prop erty of this model is the s hape of the accretion rate vs. time curve (see e.g. Fig. 4 
middle panel in iD'Angelo &: SpruitI l20ld ) . The outburst starts with a steep rise of the accretion 
rate which is followed by a slow plateau-like fading with quasi-periodic oscillations on top of the 
slow decline. This is very similar to what is observed in the optical light curve of EXLup, which 
is thought to be closely related to the accretion rate (i.e. hot spot on the stellar surface). Not 
only the shape of the light curve, but also the predicted time-scales by these simulations are in 
qualitative agreemen t with the observations. For the paramete rs which are the clo s est to that of 
EXLup (see Fig. 5 in lD'Angelo &: Spruitll20ld ) the simulations of lD'Angelo k. SpruitI (|20ld ) predict 
that the duration of the outbursts is less than the viscous time-scale at the co-rotation radius, while 
the time between two subsequent outb urst is about ten ti r nes th e viscous time-scale. To compare 
the observations to the simulation of ID'Angelo Spruit (12010 ) we estimated the viscous time- 
scale at the co-rotation radius for EXLup. Sipos et al.l ( 20od ) me asured v ■ sin(i) of EXLup t o 
be 4.4±2km/s. Using a stellar mass and radius of 0.6 1.6 R© (jGras- Velazquez &: Rayll2005l ). 
respectively, and assuming 20° inclinati on and a = 0.1 the resul ting viscous time-scale is about 



three years. According to the models of lD'Angelo Sz SpruitI (j20ld ) the duty cycle of the outbursts 
would be about thirty years while a single outburst would last less than three years. Although 
these numbers are close to the observed values, the agreement could only be considered qualitative, 
due to the large unceratainty in the value of a and the unknown inclination angle. Despite of these 



- 20 - 



uncertainties, the rema rkable similarity of the ob served light curve of EX Lup and the accretion 
rate curve predicted bv lD'Angelo &: Spruit! (j20ld ) suggests that a further quantitative comparison 
between this model and observations of EXors should be done. 



5. Summary and conclusions 

Several studies tried to find special characteristics of EXors, apart fr om the optical outbursts 



which would make them distinguishable from 'normal' TTauri stars (e.g.. lHerbie]l2007l . ISipos et al 



20091 ). but none were found so far. Thus, it was concl uded that probably most , if not all T Tauri stars 



go through an eruptive phase during their evolution (jVorobvov &: Basull2006l ) . Despite the fact that 
studies of these objects provide information about young low-mass stars in general, little is known 
about this phase of stellar evolution. In this paper we studied the disk and dust properties around 
EX Lup in its most recent outburst, in order to learn about dust processing and put constraints on 
the outburst mechanism. 

Our conclusions can be summarized as follows: 
• The accretion rate of EX Lup increased to about 2.2 • lO~^M0/yr on 2 1 Apr 2008 . Th is 



is about three orders of magnitude higher than in the quiescent phase (jSipos et al.ll2009l ). 
indicating that the outburst is caused by increased accretion onto the star. 

Modeling of the SED and mid-infrared visibilities suggests that more than 90% of the total 
accretion luminosity is radiated away from a hot component which dominates the SED at 
optical wavelengths. This component can either be a hot spot on the stellar surface or a 
narrow ring of hot gas at the inner edge of the dust disk. Our observations exclude the 
presence of an extended (Rout ~ 0.3 AU) optically thick accretion disk during the outburst. 

Calculation of the viscous time-scale in the disk shows that all material accreted in the 2008 
outburst of EX Lup should have been located within 0.1 AU from the central star. This is 
well within the inner hole of the dust disk (inner radius: 0.3 AU), suggesting that the inner 
dust-free hole is filled with gas, optically thin in the continuum. 

We compared accretion rates and time-scales as well as the SED of EX Lup and FUors. We 
concluded that, although the triggering mechanism of the outburst of EX Lup is certainly 
related to accretion, it is unlikely that it is the same thermal instability that is thought to 
operate in FUor systems. 

The post-outburst Spitzer IRS spectra of EX Lup revealed the presence of several new crys- 
talline features (19.5 ^m, 24 /im, 28 //m and tentatively at 33 /im), which are associated with 
crystalline forsterite. The detection of these forsterite features supports the conclusion of 



Abraham et al.l ( 20091 ) that crystals have formed in the disk atmosphere of EX Lup during 



the outburst. 
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On the basis of the position of the forsterite bands (especiahy the 16 /im band) we conclude 
that the lattice of the new crystals is very iron-poor. On the other hand, the crystals should 
be in thermal contact with iron, otherwise their temperature would be too low. This suggests 
that the parent amorphous grains had to be also in thermal contact with iron either in the 
silicate network or as metallic inclusions. 

Apart from forsterite, no other crystalline silicate species were identified in the mid-infrared 
spectra of EXLup. In order to prevent the formation of crystalline enstatite the amorphous 
grains should be small (<1 /im) and porous. The analysis of the quiescent phase Spitzer IRS 
spectrum of EX Lu p indeed showed that the disk atmosphere contains small amorphous grains 



(jSipos et al.ll2009l ). 



The MIDI correlated spectrum of EX Lup shows very high crystallinity, while in the uncor- 
related spectrum no crystalline features are seen. This indicates that crystals are highly 
concentrated towards the center, which is reproduced in our models as crystal formation 
occurs only within 1.1 AU from the central star. 

The shape of the uncorrelated (outer disk) spectrum shows signatures of large grains (several 
^m) in the disk atmosphere, which were not present before. This can only be explained by 
active vertical mixing, which stirred up these grains from deeper layers. It also implies that 
the turbulence in the quiescent phase is probably weaker than in the outburst, to allow these 
grains to settle below the disk atmosphere. 

Our model predicts far weaker forsterite bands longwards of 20 //m in the post-outburst 
spectra, than observed in the Spitzer spectra, indicating the presence of forsterite beyond 
1.1 AU (i.e. lower temperatures). Crystals could not form in situ outwards of 1.1 AU in the 
disk atmosphere, otherwise their feature would be visible in the MIDI uncorrelated spectrum. 
We concluded, therefore, that they were transported outwards from within 1.1 AU by some 
wind-driven transport mechanism. 

Our vertical mixing model could not reproduce the decrease in the strength of the forsterite 
bands in the 10 fim silicate complex between October 2008 and April 2009. An alternative 
explanation could be if a significant fraction of the crystal population in the innermost 1 . 1 AU 
was transported to larger radii. 
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the paper. We thank the ESO staff for executing the observations in service mode. The research of 
A. K. is supported by the Netherlands Organization for Scientific Research. The research of Zs. R. 
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A. Modeling of the outburst 



In or der to obtain the temperature structure of the disk, we used 2D radiative transfer code 
RADMC (jPullemond &: Dominikll2004al) T his code was extensively t ested against other radiative 



trans fer codes at low (jPascucci et alj 120041 ) and high optical depths (|Pinte et alj l2009l : iMin et al 



20091 ). Originally, only passive heating of the dust grains by the central star was considered in 
RADMC. For this study we also include viscous heating by an accretion disk and also the emission 
from the hot spot on the stellar surface. Viscous heating was implemented as a one grid cell layer 
in the disk midplane, which radiates as a blackbody at the local effective temperature, assuming 
it is optically thick in the vertical direction. T he effective tempera.ture of the accretion disk at a 
given radius (Tefj(r)) is calculated according to IShakura k. SunvaevI (|l973l ): 



SGM^M 



1-1^ 

r 



0.5' 



(Al) 



Here, G is the gravitational constant, M is the mass accretion rate, while and i?^ are the mass 
and radius of the star, respectively. 

We also consider the radiation of the hot spot where the funnel flow of the accreting gas 
is thought to reach the stellar surface. Assuming that the spot radiates as a blackbody, the 
temperature of the spot is given by 



GM^M 



(A2) 



Here, is the effective temperature of the star and / is the fraction of the stellar surface area 
covered by the hot spot. The quantity rtr is the truncation radius from which the accreting gas is 
thought to m ove with free-fall veloci ty towards the surface of the star. For rtr we assumed rtr = 
according to iGullbring et al.l (jl998l ). The surface area of the spot (/ = 0.5) was determined by 
fitting the o utburst SED at opti c al wa velengths. This value is higher than usually found in T Tauri 
stars (0.01, ICalvet &: Gullbring] (jl998l )). It was, however, also note d that the covering fracti on of 
the hot spot can be as high as > 0.2 in strongly accreting systems (jCalvet fc Gullbrineill998l ). To 
account for the change of the covering fraction with the accretion rate we use the following (ad 
hoc) relation for M < 2.2 • lO-'^Mo/yr: 



/ = 0.5 



M 



0.5 



2.2 • lO-^Me/yr 



(A3) 



for M > 2.2 • W-'^Mq/yi we fixed the value of / and used / = 0.5. 

After we obtained the te mperature structure of the disk at any instant, the procedure in 
DuUemond &: DominikI (|2004bl ) was followed to investigate the turbulent mixing of solids in the 
disk time-dependently. The mixing of the solid is treated as a turbulent diffusion process and we 
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also include dust settling. The conservation equation of dust grains we solve is 



df{m,z) _ d_ 
dt dz 




+ -K- {f{m,z)vsett{m,z)) = 



(A4) 



Here, f{m, z)dmdz represents the number of dust grains per square centimeter of the disk, D{m, z) 
is the diffusion coefficient and Vsett {m, z) is the settling speed. For the diffusion coefficient of the 
gas we used the standard a prescription 



where Cg is the local sound speed and Hp is the pressure scale height of the gas. The diffusion 
coefficient of the dust {D) depends on how strongly the grains are coupled to the gas (i.e. the size 
of the grains), which is described by the Schmidt-number {Sc) as D = Dq/Sc. The Schmidt number 
is unity if the dust particles are well coupled to the gas, while for larger grains sizes the decoupling 
is described by larger Schmidt numbers. During the simulation we assumed that Sc= 1. The dust 
particles, which contribute to the 10 micron emission feature of EXLup, are mostly sub-micron 
sized for which the strong coupling is a good approximation. During the calculation we assume 
that the disk is vertically isothermal and we calculate the local sound speed from the temperature 
at the disk midplane. The sound speed scales with the square root of the temperature, therefore the 
error we made with the assumption of vertically isothermal disk is in the order of (T'atm/^mid)'^'^- 
At 1 AU this factor is less than 2, which is basically negligible compared to the uncertainty in the 
value of a. 
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Table 1. Log of observations of EXLup in our campaign, supplemented with a pre-outburst 
Spitzer IRS spectrum from the Spitzer Data Archive. For photometric observations Fj, is the 
derived photometric flux density for a given filter/wavelength, while Fj^^corr is the corrected value 

for the emission lines (optical photometries only). 



Inst ruuKnil 


Dato 


A [//m] 


F„ [mJy] 


Fw ,-,„■,■ fiiiJyl 


SPITZER IRS 


18 Mar 2005 


5.5-38.0 


- 


- 


WFI 


20 Apr 2008 


0.36 (U) 


159±1.46 


159±1.46 


GROND 


20 Apr 2008 


0.45 (g) 


450±62 


293±62 


WFI 


20 Apr 2008 


0.46 (B) 


486±5 


319±5 


WFI 


20 Apr 2008 


0.54 (V) 


554±5 


419±5 


GROND 


20 Apr 2008 


0.61 (r) 


675±93 


557±93 


GROND 


20 Apr 2008 


0.77 (i) 


825±114 


754±114 


GROND 


20 Apr 2008 


0.89 (z) 


1025±141 


721±141 


SOFI 


20 Apr 2008 


1.21 


462±13 


- 


GROND 


20 Apr 2008 


1.25 (J) 


400±9 


- 


GROND 


20 Apr 2008 


1.65 (H) 


529±13 


- 


SOFI 


20 Apr 2008 


1.71 


586±16 




GROND 


20 Apr 2008 


2.15 (K) 


670±31 




SOFI 


20 Apr 2008 


2.19 


643±18 




SPITZER IRS 


21 Apr 2008 


5.5-38.0 






MIPS 


21 Apr 2008 


71.42 


3130±220 




LADOGA 


21 Apr 2008 


870 


41±10 




VLTI MIDI 


22 Jun 2008 


8.0-13.0 






VLTI MIDI 


16 Jul 2008 


8.0-13.0 






VLT VISIR 


24 Jul 2008 


8.0-13.0 






VLT VISIR 


25 Jul 2008 


8.0-13.0 






VLT VISIR 


27 Jul 2008 


8.0-13.0 






VLT VISIR 


28 Aug 2008 


8.0-13.0 






VLT VISIR 


29 Aug 2008 


8.0-13.0 






SPITZER IRS 


10 Oct 2008 


5.5-38.0 






SPITZER IRS 


7 Apr 2009 


5.5 38.0 
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Table 2. Parameters of the quiescent disk model from lSipos et alj (j2009l ). Parameters in italics 
were adopted from the literature and kept fixed, while the rest of the parameters were derived 

from fitting the SED. 



Parameter Fitted value 



Stellar temperature 3800 K 

Stellar radius 1.6 

Stellar mass 0.6 

Visual extinction Omag 

Distance 155pc 

Inner radius (dust disk) 0.3 AU 

Outer radius (dust disk) 150 AU 

Scale height (at the outer radius) 18 AU 

Flare index 0.09 

Power exponent (surface density) -1.0 

Total disk mass (gas + dust) 0.025 Mq 

Inclination angle 20 ° 
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Fig. 1. — Optical light curve of EXLup. The red line in the middle of the points shows the spline 
fit used for our 2D RT simulation. The green vertical lines mark the epochs when mid-infrared 
spectra were taken. 
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Fig. 2. — SED of EXLup in the 2008 outburst. Data shown with red are measured within 24 
hours, 21 Apr 2008. The blue hne shows the quiescent phase disk model, while the black dashed 
line shows the emission of the MO (Tefr=3800K) star in the quiescent phase. 
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Fig. 3. — Continuum-subtracted near-infrared spectrum of EXLup around the Br 7 line in the 
quiescent (black solid line) and in the outburst phase (red solid line) . The luminosity of the Br 7 
line strongly increased in the outburst compared to the quiescent phase indicating the increase of 
the accretion rate. 
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Fig. 4. — Spitzer IRS spectra of EX Lup. The red line shows the fitted continuum emission. 
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Fig. 5. — Continuum subtracted Spitzer IRS spectra of EXLup. The 5.5-13.5 micron (Left) and 
13.5-35 fiui (Right) wavelength intervals are shown separately with different scales in order to show 
also the weaker silicate band longwards of 18 /im. 




Fig. 6. — VLTI/MIDI spectra of EXLup using the U1U4 configuration with 121m projected 
baseline length. Left: Correlated spectrum (inner disk) o f EX Lup (black s o lid lin e) resembles that 
of the inner disk spectrum of HD142527 (red solid line, Ivan Boekel et al.l (|2004l )) indicating very 
high crystallinity. Right: The outer disk spectrum of EXLup is very similar to the Spitzer IRS 
spectrum of FU Ori, for which it was shown that large grains (several //m in size) are needed to 
explain the shape of the 10 //m feature complex. 
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Fig. 7. — Left: SED of EXLup in the 2008 outburst. Data shown with red are measured within 24 
hours (on 20 and 21 Apr 2008). The black dashed hue shows the MO stehar emission. Right: In 
the 1.5-20 yum wavelength interval the A-NAD models fit the observed data the best, while outside 
of this wavelength range the quality of the fit is very similar in all three model families. 
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Fig. 8. — Comparison of the observed 10 micron silicate features with om' model predictions. The 
observed Spitzer IRS spectra shown with red, while the model predictions are drawn with black. 
Each column in the figure shows the observations and result of the modeling at an epoch shown at 
the top of the columns. The different rows indicate different model series; (a) Model A-AD-aO.Ol, 
(b) Model A-AD-aO.l, (c) Model A-AD-al.O, (d) Model B-AD. 




Fig. 9. — Same as Fig. [HI but for the 13-35 /xm wavelength interval. 
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Fig. 10. — Comparison of MIDI observations (red solid lines) to our models A-AD (black solid 
lines) and A-NAD (blue solid lines). Top: Correlated (inner disk) spectrum, Middle: Uncorrelated 
(outer disk) spectrum, Bottom: Spectrally resolved visibilities. The different columns show models 
with different values of a used for the mixing simulation. The correlated spectra are reproduced 
the best with a = 0.01. The difference between our models and the observed outer disk spectra 
can be explained with the presence of large (few fim) grains in the disk atmosphere, which are not 
included in our models. Although the shape of the silicate feature is only marginally sensitive to 
the presence of an optically thick accretion disk, the visibilities are more consistent with models 
without an accretion disk. 
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a = 0,01 a = 0.1 CI = 1 .0 
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Fig. 11. — Comparison of the VISIR observations (red solid lines) to our models A- AD (black solid 
lines). Our models reproduce the strength of the silicate feature very well on 25 Jul 2008, while on 
29 Aug 2008 the model underestimates the flux. The shape of the feature is reproduced the best 
on by models with a = 0.1, however, in Aug 2008 the shape of the spectrum is somewhat uncertain 
around 9.8 ^m due to the atmospheric ozone band. 



